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Abstract 

We discuss the use of self-seeding schemes with wake monochromators to produce 
TW power, fully coherent pulses for applications at the dedicated bio-imaging 
bealine at the European X-ray FEL, a concept for an upgrade of the facility beyond 
the baseline previously proposed by the authors. We exploit the asymmetric and 
symmetric Bragg and Laue reflections (sigma polarization) in diamond crystal. 
Optimization of the bio-imaging beamline is performed with extensive start-to- 
end simulations, which also take into account effects such as the spatio-temporal 
coupling caused by the wake monochromator. The spatial shift is maximal in the 
range for small Bragg angles. A geometry with Bragg angles close to n/2 would 
be a more advantageous option from this viewpoint, albeit with decrease of the 
spectral tunability. We show that it will be possible to cover the photon energy 
range from 3 keV to 13 keV by using four different planes of the same crystal with 
one rotational degree of freedom. 



1 Introduction 

One of the highest priority for experiments at any advanced XFEL facility 
is to establish a beamline for studying biological objects at the mesoscale, 
including large macromolecules, macromolecular complexes, and cell or- 
ganelles. This requires 2-6 keV photon energy range and TW peak-power 
pulses [Hl]-|@]. However, higher photon energies (up to 13 keV) are needed 
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to reach the K-edges of commonly used elements, such as Se, for anomalous 
experimental phasing. Studies at intermediate resolutions need access to 
the water window at 0.5 keV. The pulse duration should be adjustable from 

2 fs to 10 fs. 

A basic concept and design of an undulator system for a dedicated bio- 
imaging beamline at the European XFEL was proposed in ||6j], and opti- 
mized in yZy. All the requirements in terms of photon beam characteristics 
can be satisfied by the use of a very efficient combination of self-seeding, 
fresh bunch and undulator tapering techniques laHPia]- In particular, a com- 
bination of self-seeding and undulator tapering techniques would allow to 
meet the design TW output power. The bio-imaging beamline would be 
equipped with two different self-seeding setups. For soft X-ray self-seeding, 
the monochromator consists of a grating ySy. Starting around the energy of 

3 keV it is possible to use a single crystal (wake) monochromator instead of 
a grating fl. 

In [@] we demonstrated that it will be possible to cover the photon en- 
ergy range between 3 keV and 13 keV using the C(lll), C(220) and C(004) 
symmetric Bragg reflections. In this scenario, three different crystals would 
enable self-seeding for the different spectral range. In particular, we pro- 
posed to exploit the C(lll) symmetric Bragg reflection in the photon energy 
range between 3 keV and 5 keV. 

While developing a design for the bio-imaging beamline, the authors first 
priority was to have it satisfying all requirements. Having achieved that 
goal, the next step is to optimize the design, making it as simple as possible. 
The design presented here aims for experimental simplification and cost 
reduction of the self-seeding setups. In order to improve the original design, 
here we propose to exploit asymmetric Bragg and Laue C(lll), C(113) and 
C(333) reflections together with the symmetric Bragg reflection C(004). The 
novel design of the self-seeding setup combines a wide photon energy range 
with a much needed experimental simplicity. Only one diamond crystal is 
needed, and only one rotational degree of freedom is required. 

While, in this article, we consider applications for the bio-imaging beamline 
in particular, the present study can also be applied to other beamlines, for 
example the SASE1-SASE2 lines at the European XFEL as well. 



2 Dynamical diffraction theory and asymmetric-cut crystals 

We begin our analysis of asymmetric-cut crystals by specifying the scattering 
geometry under study. With reference to Fig. 1, and following the notation 
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Fig. 1. Sketch of (a) Bragg and (b) Laue scattering geometry. 

in we identify the crystal surface with the unit vector normal to the 
crystal surface n, directed inside the crystal. The unit vector s indicates the 
trace of the Bragg planes. The direction of the incident beam is specified by 
the unit vector s , while that of the diffracted beam is given by Always 
following llzol], we call ip n the angle between n and s, xjj the angle between 
n and s , and i/^ the angle between n and s. The signs of these angles are 
univocally fixed by requiring that the angle between s and s is positive (and 
equal to the Bragg angle 6b), while the angle between s and % is negative 
(and equal to Fig. 1 shows two generic examples for both Bragg and 
Laue geometries. 

Let us consider an electromagnetic plane wave in the x-ray frequency range 
incident on an infinite, perfect crystal. Within the kinematical approxima- 
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tion, according to the Bragg law, constructive interference of waves scattered 
from the crystal occurs if the angles between the reflecting lattice planes and 
both incident and reflected beams are equal to Bragg angle Ob, Fig. 1. The 
wavelength A and the Bragg angle are related by the well known equation 

A = 2dsin(0 B ), (1) 

where d is the distance between Bragg planes, and from now on we as- 
sume reflection into the first order. This equation shows that for a given 
wavelength of the x-ray beam diffraction is possible only at certain angles 
determined by d. It is important to remember the following geometrical 
relationships: 

• In the kinematical approximation, the Bragg reflection is a mirror reflec- 
tion i.e. the angle between the incident X-ray beam and the normal to the 
Bragg planes is equal to that between the normal and the reflected beam. 

• The incident beam and the forward diffracted (i.e. the transmitted) beam 
have the same direction. 

Diffraction from an asymmetric-cut crystal is equivalent to diffraction from 



a blazed grating. This equivalence has already been noted in [1211] . The Bragg- 
reflected light from an asymmetric-cut crystal behaves as the light diffracted 
from a grating with period D-d/ cos(ip n ) and is dispersed in accordance 
with the grating equation 

A = D[sin(^ ) - sin(i^)] . (2) 

When the Bragg condition Eq. (1) and the mirror reflection condition are 
simultaneously met one obtains 

ipo = ^n + B (3) 

and 

iph = xp n - B (4) 

By direct substitution in Eq. (2) and using simple trigonometric relations, it 
is easy to check that Eq. (3) and Eq. (4) satisfy the grating equation. 

We now turn our attention beyond the kinematical approximation to the 
dynamical theory of diffraction by a crystal. With reference to Fig. 1, let 
us define the angle between the crystal surface and the input beam as 
&i - Similarly, we define the diffracted angle &d = n/2- iph, with 

6d a function of the frequency, according to Eq. (2). We also define the angle 
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between the Bragg planes and the crystal surface as a - n/2 - ip n . All these 
new quantities follow the previously defined convention as concerns their 
signs. Note that Eq. (3) and Eq. (4) can be rewritten in terms of incident and 
diffracted angles as 0/ = a - 6 B and Q d = a + 6 B . 

It is useful to describe the modification of the incident beam by means of a 
transfer function. The reflectivity curve - the reflectance - can be expressed 
in the frame of dynamical theory as 



R(di, 9 d , co) = R[Aco + co B A9i cot(0 B ), A9 d ] , (5) 

where Aco = (co-co B ),A6i = z -0 B +a,and A6 d = 6 d -0 B - oc are, respectively, 
the deviations of frequency, incident angle and diffracted (i.e. reflected) 
angle from the resonance (Bragg) frequency and angles, respectively. The 
frequency co B and the angle 6 B are given by the Bragg law co B sin(0 B ) = nc/d. 

Consider now a perfectly collimated, white beam incident on the crystal. 
Within the kinematic approximation, the transfer function as a function 
of the first argument is a Dirac <5-f unction, 5[Aco + co B A6j cot(0 B )] which is 
another representation of the Bragg law in differential form: 



^- = Acot(0 B ). (6) 

In addition to this, within the kinematic approximation, reflection from a 
lattice plane is always a mirror reflection, and 



Moving to the framework of the dynamical theory, and in contrast to what 
has just been said, the reflectivity width is finite. This means that there is 
a reflected beam even when the incident angle and the wavelength of the 
incoming beam are not exactly related by the Bragg equation. It is interesting 
to note that incident beam and transmitted beam continue to have the 
same direction also in the framework of dynamical theory. However, in 
the framework of the dynamical theory and in the case of asymmetric-cut 
crystal, the reflection from lattice planes are not a mirror reflection anymore: 
the grating equation Eq. (2), holds, and the reflection from the crystal surface 
is always a grating reflection. 

The diffraction from the asymmetric-cut crystal in Bragg and Laue geometry 
can be described an asymmetry parameter, b, defined by 
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n ■ s h cos(xp h ) ' 

In literature on X-ray crystal diffraction, the asymmetry parameter is some- 
times defined as inverse of b, and indicated with the letter y, for example in 
Here we follow the same convention used in ||21fl. Note that \b\ = o in /o oui 
is the ratio of the widths of incident and diffracted beams. From Eq. (2) fol- 
lows that 



As has been pointed out elsewhere, this is a consequence of Liouville the- 
orem P21Q . In fact, a quasi-monochromatic, collimated beam of finite width 
that is reduced in width must be dispersed in angle to preserve the space- 
angle phase volume. 

For the Bragg geometry, where the diffracted light exists on the same side 
of the surface as the incident beam, we have b < 0, whereas b > describes 
Laue geometry. In particular, b = -1 corresponds the symmetric Bragg 
diffraction, and b = 1 corresponds to the symmetric Laue diffraction. 

The dispersion of the diffracted light from the asymmetric-cut crystal fol- 
lows directly from the grating equation, yielding 



A— £ = (l + &)tan(0 B ). (10) 

uA 



As expected, there is no dispersion for symmetric Bragg geometry (b = -1) 
because there is no grating structure in that case. 

In the following, we will consider one single diamond crystal, as the 100 /im 
thick crystal proposed in [14], and currently used into the LCLS self-seeding 
setup ||22|]. We will define a cartesian reference system {x, y,z) linked with 
the crystal. The direction z corresponds to the direction identified by the 
Miller indexes 0, 0, 1, which also coincides with the direction of n, while x 
and y are specified as in Fig. 2, so that the direction [110] corresponds to the 
x direction. 

The crystal can rotate freely around the y axis (pitch angle) as indicated in 
the figure. In this way, we can exploit several symmetric and asymmetric 
reflection: we will be able to cover the entire energy range between 3 keV 
and 13 keV by changing the pitch angle of the crystal in Fig. 2. In the low 
energy range between 3 keV and 5 keV we will use the C(lll) asymmetric 
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Fig. 2. Drawing of the orientation of the diamond crystal considered in this article. 



Table 1 

Useful Diamond Reflections 



hkl 


Min. Energy (keV) 


AE (meV) 


111 


3.01034 


192.0 


311 


5.76401 


56.0 


400 


6.95161 


60.6 


333 


9.03035 


27.3 


444 


12.0404 


25.0 



reflection (in Bragg or Laue geometry, depending on the energy). At higher 
energies between 5 keV and 7 keV we propose to go to the C(113) asymmetric 
reflection in Bragg geometry. In the range between 7 keV and 9 keV we can 
use instead the C(004) symmetric reflection, in Bragg geometry. Finally, for 
energies larger than 9 keV wepropose the use of the C(333) asymmetric 
reflection, in Bragg geometryli]. The list of reflections considered in this 
article is summarized in Table 1, extracted from [1230 . 

2 It should be remarked that above 12 keV the C(444) asymmetric Bragg reflection 
turns out to be a valid alternative, though results are not explicitly analyzed in 
this article. The Bragg angle is, in this case, very near to n/2, 6b = 84.4 deg, while 
ip n = 35.3 deg. The relative Darwin width amount to about 2 • 10" 6 , with a very high 
transmissivity (about 97%). Our calculations show that one can reach pulse power 
in excess of 600 GW at 12 keV with a FWHM relative bandwidth of about 5 • 10" 5 . 
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Fig. 3. Scattering geometry, modulus and phase of the transmittance for the C(lll) 
asymmetric Bragg reflection from the 100 fim thick perfect diamond crystal in Fig. 

2 at 3.5 keV. 

For self-seeding implementation, we are interested in the forward diffracted, 
i.e. in the transmitted beam for each of these reflections. From this viewpoint, 
the crystal can be characterized as a filter with given complex transmissivity. 
We will now consider the three transmissivity separately, showing their 
amplitude and phase, together with the crystal geometry for fixed energy 
points. 

In Fig. 3 we show the amplitude and the phase for the C(lll) asymmetric 
Bragg reflection at 3.5 keV. The following Fig. 4 shows the same reflection 
at a higher energy, 5 keV. In this case, one has Laue geometry. In Fig. 5 we 
plot the same quantities for the C(113) asymmetric Bragg geometry at 7 keV. 
In Fig. 6 we consider the symmetric Bragg reflection C(004) at 8 keV and, 
finally, in Fig. 7 we plot amplitude and phase for the asymmetric Bragg 
reflection C(333) at 9 keV. By inspecting the plots, the reader can easily see 
the difference between Bragg and Laue geometry, as well as that between 
symmetric and asymmetric reflections. 

Amplitude and phase in the plots are related by Kramers-Kroning re- 
lations. Let us regard our crystal as a filter with transmission T(co) = 
\T(co)\ exp[zO(o;)]. According to Titchmarsch theorem (see [24] and refer- 
ences therein for a recent review on the subject) causalityLj and square- 

3 Causality simply requires that the filter can respond to a physical input after the 
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Fig. 4. Scattering geometry, modulus and phase of the transmittance for the asym- 
metric C(lll) Laue reflection from the 100 fim thick perfect diamond crystal in Fig. 
2 at 5 keV. 

integrability of the inverse Fourier transform of T(co) - \T(co)\ exp[iQ>(co)], 
which will be indicated with T{t), is equivalent to the existence of an an- 
alytic continuation of T(co) to Q = co + ico' on the upper complex Q-plane 
(i.e. for co' > 0). The same theorem also shows that the two previous state- 
ments are equivalent to the fact that real and imaginary part of T(co) are 
connected by Hilbert transformation. Since T{t) must be real (thus imply- 
ing that T*(co) - T(-co)), from the Hilbert transformation follows the well- 
known Kramers-Kroninig relations ijz^El, linking real and imaginary part 
of T(co). 

A similar reasoning can be done for the modulus \T(co)\ and the phase Q(co), 
see [27]. In fact, one can write 



ln[TM] = M\T(G))\] + iO(co) , (11) 

that play a similar role to real and imaginary part of the refractive index of 
a given medium. 

Note that T{co) = T(-co) implies that \T(co)\ = \T(-co)\ and that <D(<y) = 
-0(-co). Therefore, using Eq. (11) one also has that ln[T(a;)]* = ln[T(-o;)]. 

time of that input and never before. 

4 Here we are tacitly assuming, as required by Titchmarsh theorem, that T(co) is a 
square integrable function. 
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C{113)at7keV 
Bragg geometry 



Abs[T] 




0,1778225 0,1778250 0,1778275 0,1778300 
k [nm] 



Fig. 5. Scattering geometry, modulus and phase of the transmittance for the asym- 
metric C(113) Bragg reflection from the 100 fim thick perfect diamond crystal in 
Fig. 2 at 7 keV. 

Application of Titchmarsh theorem shows that the analyticity of ln[|T(Q)|] 
on the upper complex Q-plane implies that O(co) can be obtained from |T(o;)|. 
As is well known however, in applying such procedure one tacitly assumes 
that ln[T(Q)] is analytical on the upper complex Q-plane. While causality 
implies this fact for T(Q), it does not imply it automatically for ln[|T(Q)|]. 
In fact, such function is singular where T(Q) = 0. If T(Q) has zeros on the 
upper complex plane, these zeros would contribute adding extra terms to 
the total phase. It can be shown that T(Q) does not have zeros on the upper 
complex plane, and that Titchmarsch theorem actually applies. 

The simplest check of the correctness of the calculations by means of dy- 
namical theory consists, therefore, in verifying that amplitude and phase of 
the transmissivity are actually related by Kramer s-Kroning relations. 



3 Spatiotemporal shift and its influence on input coupling factor 



A wake monochromator introduces spatiotemporal deformations of the 
seeded X-ray pulse, which can be problematic for seeding. The spatiotem- 
poral coupling in the electric field relevant to self-seeding schemes with 
crystal monochromators has been analyzed in the framework of classical 
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I [nm] 

Fig. 6. Scattering geometry, modulus and phase of the transmittance for the sym- 
metric C(004) Bragg reflection from the 100 fim thick perfect diamond crystal in 
Fig. 2 at 8 keV. 

dynamical theory of X-ray diffraction [ 123Q . 

This analysis shows that a crystal in Bragg or Laue geometry transforms 
the incident electric field E(x, t) in the {x, f } domain into E(x - at, £), where 
a - const. The physical meaning of this distortion is that the beam spot 
size is independent of time, but the beam central position changes as the 
pulse evolves in time. First we will show in a simple manner that, based on 
the only use of the Bragg law, we can directly arrive to an explanation of 
spatiotemporal coupling phenomena in the dynamical theory of diffraction 

Q. 

The transmissivity curve - the transmittance - in asymmetric and symmetric 
Bragg or Laue geometry can be expressed in the framework of dynamical 
theory with the help of Eq. (5). Transmitted (i.e. forward diffracted) beam 
and incident beam have the same direction. Therefore we set (9 Z = 6^ Using 
the symbol T for transmittance, Eq. (5) yields 

T(Aco, A9i) = T(Aco + co B A9 cot(0 B )) , (12) 

where Aco - (co - co B ) and AO - (6 l ■ - 6b) are the deviations of frequency 
and incident angle of the incoming beam from Bragg frequency and Bragg 
angle respectively. The frequency co B and the angle 9 B are obviously related 
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Fig. 7. Scattering geometry, modulus and phase of the transmittance for the asym- 
metric C(333) Bragg reflection from the 100 fim thick perfect diamond crystal in 
Fig. 2 at 9 keV. 

by the Bragg law: co B sin(0 B ) = nc/d. 

It should be realized that the crystal does not introduce an angular disper- 
sion similar to the grating. However, a more detailed analysis based on the 
expression for the transmissivity, Eq. (12), shows that a less well-known 
spatiotemporal coupling exists. Let us discuss this fact. It is evident by in- 
spection of Eq. (12), that the transmissivity is invariant under angle and 
frequency transformations obeying 

Aco + co B AG cot(0 B ) = const . (13) 

This corresponds to the coupling in the Fourier domain. In general, one 
would indeed expect the transformation to be symmetric in both the Fourier 
domain {k x ,co} (with k x - co B A6/c) and in the space-time domain {x, £}, 
due to the symmetry of the transfer function. However, it is reasonable to 
expect the influence of a nonsymmetric input beam distribution. The field 
transformation for the XFEL pulse after the crystal in the \x, t] domain is 
given by 



where FT indicates a Fourier transform from the frequency domain [k x/ co) 



Eout(*, t) = FT{T[Ag), (k x ) in , (k x ) ou ^E^(Aa) f k x )} , 



(14) 
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to the space-time domain {x, £}. 



In the self-seeding case, the incoming XFEL beam is well collimated, mean- 
ing that its angular spread is a few times smaller than the angular width of 
the transfer function. Yet, the spectral bandwidth of the incoming beam is 
much wider than the bandwidth of the transfer function. Under the limit 
of a wide spectral bandwidth for the incoming beam (with respect to the 
bandwidth of the transfer function), and applying the shift theorem twice 
we obtain 



E Q ut(x,t) = t](t) J dk x exp(i cot 6 B k x ct)exp(-ik x x)E(co B ,k x ) 

= t](t)a(x-c cot 6 B t), (15) 



where 



7](t) =2^J dcoexp(-icot)T(co) (16) 

is the inverse temporal Fourier transform of the transmittance curve (im- 
pulse response). Eq. (15) is universal i.e. valid for both symmetric and asym- 
metric Bragg and Laue geometries. 

The spatial shift given by Eq. (15) is proportional to cot(0 B ), and is maximal in 
the range for small values of 6 B . Therefore, the spatiotemporal coupling is an 
issue, and efforts are necessarily required to avoid distortions. A geometry 
with 6 B close to tz/2 would be a more advantageous option from this view 
point, albeit with decrease in the spectral tunability [23]. 



3.1 Influence of the offset on the coupling factor 



The influence of the spatiotemporal distortion on the operation of the self- 
seeding setup can be quantified by studying the input coupling factor be- 
tween the seed beam and the ground eigenmode of the FEL amplifier. When 
the undulator is sufficiently long, the output energy in the FEL pulse grows 
exponentially with the undulator length, and the eenrgy gain, G = E out /E ext , 
can be written as 



G = Aexp[z/L g ], (17) 

where A is the input coupling factor. In the linear regime the energy gain 
does not depend on the input energy E ext/ and for the case including spatial 
shift the input coupling factor A allows for a simple and convenient measure 
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shifted 



0,1512 



x[mm] S [nm] X[nm] 



Fig. 8. Comparison between the self-seeded output with (red, dotted lines) and 
without (black, solid lines) transverse offset, at different offset values. The left plots 
show the photon beam transverse profiles, the middle plots show the power, the 
right plots show the spectra. All plots refer to a single simulation shot. 

of the departure from the ideal situation. In our case this measure is simply 
defined as 



q _ ^nonideal _ ^-nonideal (18) 
^ ideal ^-ideal 

We will consider a particular example for A = 0.15 nm, corresponding to 
about 8 keV. We simulated the SASE process in the first 11 cells, followed 
by the self-seeding setup exploiting the C(004) reflection. The seed beam is 
then superimposed again onto the electron bunch (the modulation having 
been washed out by the chicane). We calculate the output power and spec- 
trum after 6 cells, which is still in the linear regime, and we scan over the 
transverse offset. 

A few examples of the output obtained with this procedure are shown in Fig. 
8, which compares the result for zero offset (black solid lines) with those for 
2 jum, 18 /im, and 34 /im (red dotted lines). All results refer to a single sim- 
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9,0x1 0~ 5 n 




z[m] 

Fig. 9. Collected energies as a function of the position inside the radiator, for 
increasing offset from fim to 40 fim with steps of 2 /im. 




0,0 -\ — i — i — i — i — i — i — i — i — i — i — i — i — i — 
5 10 15 20 25 30 

A[]um] 

Fig. 10. Measure of the dependence of the input coupling factor on the offset Q(A). 

ulation shot. The plots on the left show the transverse profile of the photon 
beam after 6 undulator cells following the self-seeding setup. The figures in 
the middle show the power, and those on the left the spectra. Qualitatively, 
the effect due to the transverse offset is easy to see by inspection. 

The energy as a function of the position into the radiator is shown in Fig. 9. 
The function Q(A) instead, is shown in Fig. 10. Note that the input coupling 
factor does not drop abruptly to zero when A is of order of the rms transverse 
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Fig. 11. Schematic setup of the bio-imaging beamline discussed in |@]. 

size of the photon beam. This is because the photon beam intensity does not 
drop off as a Gaussian function. As a result, on can qualitatively go to larger 
deviations of 6b , still obtaining a reasonable coupling factor. 

We hope that these considerations will be useful for FEL physicists in the 
design stage of wake monochromator setups. 



4 FEL studies 



In order to quantitatively investigate the effects of the spatial shift on the 
FEL amplification process, we performed numerical simulations using the 
code Genesis 129Q . Simulations are based on a statistical analysis consisting 
of 100 runs. We consider the setup in Fig. 11, with reference to |@] for details. 
In this article we are interested in the energy range where hard X-ray self 
seeding can be implemented, starting from 3.5 keV up to 13 keV. The C(lll) 
reflections (Bragg and Laue) will be used in the second chicane. C(113), 
C(004) and C(333) reflections will be used in the third chicane. Since we are 
interested in energies starting from 3.5 keV, we fix the electron energy at 
17.5 GeV. The lower energy is used in the very soft X-ray regime, see [0], 
between 0.3 keV and 0.5 keV, which is outside the interests of this article, and 
where one relies on a grating monochromator. In Fig. 12 we plot the results 
of start-to-end simulations 11301 ], while the main undulator parameters are 
reported in Table 2. Our simulations automatically include the influence of 
the spatial shift of the photon beam induced by the crystal. We will present 
results for the different reflections at given energy points. 

As is explained in detail in [7], in the range between 3 keV and 5 keV the 
first chicane is not used and is switched off. After the first 7 cells the electron 
and the photon beams are separated with the help of the second magnetic 
chicane, and the C(lll) reflection is used to monochromatize the radiation. 
The seed is amplified in the next 4 cells. After that, the electron and the 
photon beam are separated again by the third chicane, and an X-ray optical 
delay line allows for the introduction of a tunable delay of the photon beam 
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Fig. 12. Results from electron beam start-to-end simulations at the entrance of 
the undulator system of the bio-imaging beamline [301] for the hard X-ray case 
for the 17.5 GeV mode of operation. (First Row, Left) Current profile. (First Row, 
Right) Normalized emittance as a function of the position inside the electron beam. 
(Second Row, Left) Energy profile along the beam, lower curve. The effects of 
resistive wakefields along SASE1 are illustrated by the comparison with the upper 
curve, referring to the entrance of SASE1 (Second Row, Right) Electron beam energy 
spread profile, upper curve. The effects of quantum diffusion along SASE1 are 
illustrated by the comparison with the lower curve, referring to the entrance of 
SASE1. (Bottom row) Resistive wakefields in the SASE3 undulator [301] . 

with respect to the electron beam. The following 6 cells use only a part of 
the electron beam as a lasing medium. A magnetic chicane follows, which 
shifts the unspoiled part of the electron bunch on top of the of the photon 
beam. In this way, a fresh bunch technique can be implemented. Since the 
delays are tunable, the photon pulse length can be tuned. Finally, radiation 
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Table 2 

Undulator parameters 



Units 



Undulator period mm 68 
Periods per cell - 73 

Total number of cells - 40 
Intersection length m 1.1 

Photon energy keV 0.3-13 



„ 1,5x10 lJ 
Z> 




6,0x1 9 



0,355 0,356 

X[nm] 




Fig. 13. Power and spectrum at 3.5 keV before the second magnetic chicane. Grey 
lines refer to single shot realizations, the black line refers to the average over a 
hundred realizations. 

is amplified into the last 23 tapered cells to provide pulses with about 2 TW 
power. 

In partial difference with respect to what has been discussed in [7], for 
energies lower than 5 keV we used the asymmetric C(lll) reflection (Bragg 
and Laue geometry), while in [7] we exploited the C(lll) symmetric Bragg 
reflection only. Moreover, for energies larger than 5 keV we rely on a crystal 
placed after 11 cells, and we take advantage of the C(113), the C(004) and the 
C(333) reflections. In this case the fourth chicane is switched off. If tunability 
of the pulse duration is requested in this energy range, this is most easily 
achieved by providing additional delay with the fourth magnetic chicane 
installed behind the hard X-ray self-seeding setup. 



4.1 C(lll) asymmetric Bragg reflection at 3.5 keV 



As discussed above, the first chicane is switched off, so that the first part of 
the undulator effectively consists of 7 uniform cells. We begin our investiga- 
tion by simulating the SASE power and spectrum after the first part of the 
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Fig. 14. Power and spectrum at 3.5 keV after the single crystal self-seeding 
X-ray monochromator. A 100 fim thick diamond crystal ( C(lll) Bragg reflection, 
cr-polarization ) is used. Grey lines refer to single shot realizations, the black line 
refers to the average over a hundred realizations. The black arrow indicates the 
seeding region. 
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Fig. 15. Logarithmic plot of the power at 3.5 keV after the single crystal self-seeding 
X-ray monochromator. The black arrow indicates the seeding region. The region 
on the right hand side of the plot (before the FEL pulse) is nominally zero because 
of causality reasons. Differences with respect to zero give back the accuracy of our 
calculations. 

undulator, that is before the second magnetic chicane in the setup. Results 
are shown in Fig. 13. 

The second magnetic chicane is switched on, and the single-crystal X-ray 
monochromator is set into the photon beam. We take advantage of the C(l 1 1) 
Bragg reflection, a-polarization, Fig. 14. As discussed before, the modulus 
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Fig. 16. Power and spectrum at 3.5 keV after the third chicane equipped with the 
X-ray optical delay line, delaying the radiation pulse with respect to the electron 
bunch. Grey lines refer to single shot realizations, the black line refers to the average 
over a hundred realizations. 

and the phase of the crystal transmissivity are related by Kramers-Kroening 
relations. The numerical accuracy with which causality is satisfied in our 
simulations can be shown by a logarithmic plot of the FEL pulse power after 
the crystal, Fig. 15. The peak in the center is the main FEL pulse. On the left 
side one can identify the seed pulse. The large jump between the power after 
transmission (on the right side of the peak in the center) is due to causality, 
and is strictly connected with the fact that the phase in the transmission 
function is properly accounted for. For example, if a transmission without 
phase were used, one would have obtained a symmetric behavior, which 
is not causal at all. The small (~ 10" 3 ) but visible departure from the exact 
zero appearing in the right-hand side of Fig. 15 is related with the accuracy 
of our calculations. This accuracy is acceptable for most purposes. 

Following the seeding setup, the electron bunch amplifies the seed in the 
following 4 undulator cells. After that, a third chicane is used to allow for 
the installation of an x-ray optical delay line, which delays the radiation 
pulse with respect to the electron bunch. The power and spectrum of the 
radiation pulse after the optical delay line are shown in Fig. 16, where the 
effect of the optical delay is illustrated. 

Due to the presence of the optical delay, only part of the electron beam is 
used to further amplify the radiation pulse in the following 6 undulator 
cells. The electron beam part which is not exploited is fresh, and can further 
lase. In order to allow for that, after amplification, the electron beam passes 
through the final magnetic chicane, which delays the electron beam. The 
power and spectrum of the radiation pulse after the last magnetic chicane 
are shown in Fig. 17. By delaying the electron bunch, the magnetic chicane 
effectively shifts forward the photon beam with respect to the electron beam. 
Tunability of such shift allows the selection of different photon pulse length. 
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Fig. 17. Power and spectrum at 3.5 keV after the last magnetic chicane. Grey lines 
refer to single shot realizations, the black line refers to the average over a hundred 
realizations. 
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Fig. 18. Tapering law at 3.5 keV. 

The last part of the undulator is composed by 23 cells. It is partly tapered 
post-saturation, to increase the region where electrons and radiation interact 
properly to the advantage of the radiation pulse. Tapering is implemented by 
changing the K parameter of the undulator segment by segment according 
to Fig. 18. The tapering law used in this work has been implemented on an 
empirical basis. 

The use of tapering together with monochromatic radiation is particularly 
effective, since the electron beam does not experience brisk changes of the 
ponderomotive potential during the slippage process. The final output is 
presented in Fig. 19 in terms of power and spectrum. As one can see, simu- 
lations indicate an output power of about 2 TW. 

The energy of the radiation pulse and the energy variance are shown in Fig. 
20 as a function of the position along the undulator. The divergence and 
the size of the radiation pulse at the exit of the final undulator are shown, 
instead, in Fig. 21. In order to calculate the size, an average of the transverse 
intensity profiles is taken. In order to calculate the divergence, the spatial 
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Fig. 19. Final output. Power and spectrum after tapering at 3.5 keV. Grey lines 
refer to single shot realizations, the black line refers to the average over a hundred 
realizations. 




Fig. 20. Final output. Energy and energy variance of output pulses at 3.5 keV. In 
the left plot, grey lines refer to single shot realizations, the black line refers to the 
average over a hundred realizations. 

Fourier transform of the field is calculated. 



4.2 C(lll) asymmetric Laue reflection at 5.0 keV 

Operation at 5 keV is identical to the case for 3.5 keV. The only difference is 
that now the Laue C(lll) reflection is used. 

As before, the first chicane is switched off, so that the first part of the un- 
dulator effectively consists of 7 uniform cells. After the first part of the 
undulator, one has the power and spectrum shown in Fig. 22. The second 
magnetic chicane is switched on, and the single-crystal X-ray monochro- 
mator is set into the photon beam. We take advantage of the C(lll) Laue 
reflection, a-polarization, Fig. 23. 
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Fig. 21. Final output. X-ray radiation pulse energy distribution per unit surface 
and angular distribution of the X-ray pulse energy at 3.5 keV at the exit of output 
undulator. 
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Fig. 22. Power and spectrum at 5 keV before the second magnetic chicane. Grey 
lines refer to single shot realizations, the black line refers to the average over a 
hundred realizations. 

Exactly as before, following the seeding setup, the electron bunch amplifies 
the seed in the following 4 undulator cells. After that, a third chicane is used 
to allow for the installation of an x-ray optical delay line, which delays the 
radiation pulse with respect to the electron bunch. The power and spectrum 
of the radiation pulse after the optical delay line are shown in Fig. 25, where 
the effect of the optical delay is illustrated. The numerical accuracy with 
which causality is satisfied in our simulations can be shown by a logarithmic 
plot of the FEL pulse power after the crystal, Fig. 24. The peak in the center 
is the main FEL pulse. On the left side one can identify the seed pulse. One 
the right side, before the FEL pulse, one has nominally zero power. 

Due to the presence of the optical delay, only part of the electron beam is 
used to further amplify the radiation pulse in the following 6 undulator 
cells. The electron beam part which is not used is fresh, and can be used for 
further lasing. In order to do so, after amplification, the electron beam passes 
through the final magnetic chicane, which delays the electron beam. The 
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Fig. 23. Power and spectrum at 5 keV after the single crystal self-seeding 
X-ray monochromator. A 100 fim thick diamond crystal ( C(lll) Laue reflection, 
cr-polarization ) is used. Grey lines refer to single shot realizations, the black line 
refers to the average over a hundred realizations. The black arrow indicates the 
seeding region. 
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Fig. 24. Logarithmic plot of the power at 5 keV after the single crystal self-seeding 
X-ray monochromator. The black arrow indicates the seeding region. The region 
on the right hand side of the plot (before the FEL pulse) is nominally zero because 
of causality reasons. Differences with respect to zero give back the accuracy of our 
calculations. 

power and spectrum of the radiation pulse after the last magnetic chicane 
are shown in Fig. 26. By delaying the electron bunch, the magnetic chicane 
effectively shifts forward the photon beam with respect to the electron beam. 
Tunability of such shift allows the selection of different photon pulse length. 

The last part of the undulator is composed by 23 cells. It is partly tapered 
post-saturation, to increase the region where electrons and radiation interact 
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Fig. 25. Power and spectrum at 5 keV after the third chicane equipped with the 
X-ray optical delay line, delaying the radiation pulse with respect to the electron 
bunch. Grey lines refer to single shot realizations, the black line refers to the average 
over a hundred realizations. 
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Fig. 26. Power and spectrum at 5 keV after the last magnetic chicane. Grey lines 
refer to single shot realizations, the black line refers to the average over a hundred 
realizations. 
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Fig. 27. Tapering law at 5 keV. 

properly to the advantage of the radiation pulse. Tapering is implemented by 
changing the K parameter of the undulator segment by segment according 
to Fig. 27. The tapering law used in this work has been implemented on an 
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Fig. 28. Final output. Power and spectrum after tapering at 5 keV. Grey lines refer 
to single shot realizations, the black line refers to the average over a hundred 
realizations. 




Fig. 29. Final output. Energy and energy variance of output pulses at 5 keV. In 
the left plot, grey lines refer to single shot realizations, the black line refers to the 
average over a hundred realizations. 



empirical basis. 

As explained before, the use of tapering together with monochromatic ra- 
diation is particularly effective. The final output is presented in Fig. 28 in 
terms of power and spectrum. As one can see, simulations indicate an output 
power of about 2 TW. 

The energy of the radiation pulse and the energy variance are shown in Fig. 
29 as a function of the position along the undulator. The divergence and 
the size of the radiation pulse at the exit of the final undulator are shown, 
instead, in Fig. 30. 
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Fig. 30. Final output. X-ray radiation pulse energy distribution per unit surface 
and angular distribution of the X-ray pulse energy at 5 keV at the exit of output 
undulator. 
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Fig. 31. Power and spectrum at 7 keV before the third magnetic chicane. Grey lines 
refer to single shot realizations, the black line refers to the average over a hundred 
realizations. 

4.3 C(113) asymmetric Bragg reflection at 7 keV 



We now consider the energy point at 7 keV. In this case a simpler scheme with 
a single hard x-ray self-seeding setup is considered. The first two magnetic 
chicanes are switched off. Therefore, the electron beam lases in SASE mode 
along the first 11 undulator cells before passing through the single-crystal 
monochromator filter. The input power and spectrum are shown in Fig. 31. 

The effect of the filtering process is illustrated, both in the time and in the 
frequency domain, in Fig. 32. 

The numerical accuracy with which causality is satisfied in our simulations 
can be shown by a logarithmic plot of the FEL pulse power after the crystal, 
Fig. 33. The peak in the center is the main FEL pulse. On the left side one 
can identify the seed pulse. One the right side, before the FEL pulse, one 
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Fig. 32. Power and spectrum at 7 keV after the single crystal self-seeding X-ray 
monochromator. A 100 fim thick diamond crystal ( C(113) Bragg reflection, 
cr-polarization ) is used. Grey lines refer to single shot realizations, the black line 
refers to the average over a hundred realizations. The black arrow indicates the 
seeding region. 




shim] 

Fig. 33. Logarithmic plot of the power at 7 keV after the single crystal self-seeding 
X-ray monochromator. The black arrow indicates the seeding region. The region 
on the right hand side of the plot (before the FEL pulse) is nominally zero because 
of causality reasons. Differences with respect to zero give back the accuracy of our 
calculations. 

has nominally zero power. 

The seed is amplified up to saturation in the output undulator. As already 
discussed, we can use post-saturation tapering to increase the output power 
level. The tapering configuration in Fig. 34 is optimized for maximum output 
power level. 
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Fig. 34. Tapering law at 7 keV. 
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Fig. 35. Final output in the case of tapered output undulator at 7 keV. Power and 
spectrum are shown. Grey lines refer to single shot realizations, the black line refers 
to the average over a hundred realizations. 




Fig. 36. Energy and energy variance of output pulses in the case of tapered output 
undulator at 7 keV. In the left plot, grey lines refer to single shot realizations, the 
black line refers to the average over a hundred realizations. 

The output characteristics, in terms of power and spectrum, are shown in 
Fig. 35. The output power is increased of about a factor 20, allowing one 
to reach about one TW. The spectral width remains almost unvaried. The 
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Fig. 37. Final output. X-ray radiation pulse energy distribution per unit surface 
and angular distribution of the X-ray pulse energy at 7 keV at the exit of output 
undulator. 

output level has been optimized by changing the tapering law, resulting 
in Fig. 34, and by changing the electron beam transverse size along the 
undulator, as suggested in Jl9[|. Optimization was performed empirically. 
The evolution of the energy per pulse and of the energy fluctuations as a 
function of the undulator length are shown in Fig. 36. Finally, the transverse 
radiation distribution and divergence at the exit of the output undulator are 
shown in Fig. 37. 

4.4 C(004) symmetric Bragg reflection at 8 keV 

A feasibility study dealing with this energy range can be found in |@]. No 
changes are foreseen in this case. 

4.5 C(333) asymmetric Bragg reflection at 9 keV 

Finally, we consider the energy point at 9 keV, which relies on the C(333) 
asymmetric Bragg reflection. Again, the first two magnetic chicanes are 
switched off. The electron beam lases in SASE mode along the first 11 undu- 
lator cells before passing through the single-crystal monochromator filter. 
Input power and spectrum are shown in Fig. 38 

The effect of the filtering process is illustrated, both in the time and in the 
frequency domain, in Fig. 39. 

The numerical accuracy with which causality is satisfied in our simulations 
can be shown by a logarithmic plot of the FEL pulse power after the crystal, 
Fig. 40. The peak in the center is the main FEL pulse. On the left side one 
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Fig. 38. Power and spectrum at 9 keV before the third magnetic chicane. Grey lines 
refer to single shot realizations, the black line refers to the average over a hundred 
realizations. 
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Fig. 39. Power and spectrum at 9 keV after the single crystal self-seeding X-ray 
monochromator. A 100 fim thick diamond crystal ( C(333) Bragg reflection, 
cr-polarization ) is used. Grey lines refer to single shot realizations, the black line 
refers to the average over a hundred realizations. The black arrow indicates the 
seeding region. 

can identify the seed pulse. One the right side, before the FEL pulse, one 
has nominally zero power. 

The seed is amplified up to saturation in the output undulator. As already 
discussed, we can use post-saturation tapering to increase the output power 
level. The tapering configuration in Fig. 41 is optimized for maximum output 
power level. 

The output characteristics, in terms of power and spectrum, are shown in 
Fig. 42. The output power is increased of about a factor 20, allowing one 
to reach about one TW. The spectral width remains almost unvaried. The 
output level has been optimized by changing the tapering law, resulting 
in Fig. 41, and by changing the electron beam transverse size along the 
undulator, as suggested in |19l ]. Optimization was performed empirically. 
The evolution of the energy per pulse and of the energy fluctuations as a 
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Fig. 40. Logarithmic plot of the power at 9 keV after the single crystal self-seeding 
X-ray monochromator. The black arrow indicates the seeding region. The region 
on the right hand side of the plot (before the FEL pulse) is nominally zero because 
of causality reasons. Differences with respect to zero give back the accuracy of our 
calculations. 
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Fig. 41. Tapering law at 9 keV. 

function of the undulator length are shown in Fig. 43. Finally, the transverse 
radiation distribution and divergence at the exit of the output undulator are 
shown in Fig. 44. 



5 Conclusions 



Self-seeding scheme with wake monochromators are routinely used in gen- 
eration of narrow bandwidth X-ray pulses at LCLS (220 . Recently, the pho- 
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Fig. 42. Final output in the case of tapered output undulator at 9 keV. Power and 
spectrum are shown. Grey lines refer to single shot realizations, the black line refers 
to the average over a hundred realizations. 
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Fig. 43. Energy and energy variance of output pulses in the case of tapered output 
undulator at 9 keV. In the left plot, grey lines refer to single shot realizations, the 
black line refers to the average over a hundred realizations. 

ton energy range of self-seeding setup at LCLS was extended from 8 keV 
to lower photon energies down to 5.5 keV by exploiting different planes on 
the same diamond crystal [13 II ]. 

In this article we propose a study of the performance of self -seeding setups 
for the European XFEL based on wake monochromators, exploiting dif- 
ferent (symmetric and asymmetric, Bragg and Laue) reflections of a single 
diamond crystal, similar to that used into the LCLS setup. In our study we 
account the spatiotemporal shift of the seed induced by the passage through 
the crystal, and we study the dependence of the FEL input coupling factor 
on this shift. We exemplify our results using the concept for a dedicated 
bio-imaging beamline previously proposed by the authors. We conclude 
that crystals similar to the one actually used at the LCLS can be exploited to 
cover the energy range between 3 keV and 13 keV with high performance, 
allowing one to reach TW-class pulses. 
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Fig. 44. Final output. X-ray radiation pulse energy distribution per unit surface 
and angular distribution of the X-ray pulse energy at 9 keV at the exit of output 
undulator. 
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